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Q ■ The scalar field theory with higher derivatives is considered in the 

first order formalism. The field equation of the forth order describes 
scalar particles possessing two mass states. The first order relativis- 
tic wave equation in the 10-dimensional matrix form is derived. We 
find the relativistically invariant bilinear form and corresponding La- 
£NJ I grangian. The canonical energy- momentum tensor and density of the 

electromagnetic current are obtained. Dynamical and non-dynamical 
. components of the wave function are separated and the quantum- 

mechanical Hamiltonian is found. Projection operators extracting so- 
lutions of field equations for definite energy and different mass states 
of particles are obtained. The canonical quantization of scalar fields 
with two mass states is performed, and propagators are found in the 
formalism considered. 
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^ 1 Introduction 

It is known that the gravity theory based on the Einstein-Hilbert action is 
non-renormalizable in four dimensions pLj. By including curvature squared 
terms in the action [2], the theory becomes renormalizable but the higher 
derivative (HD) theory. There are also HD fields in generalized electrody- 
namics [3j. HD theories allow us to improve renormalization of the theory 
and to avoid ultraviolet divergences pE]. However, it was discovered soon that 
there are some difficulties with negative norm (ghosts), and unitarity in HD 
theories [5], [6]. Much attention, therefore, was made to study the simplest 
version of the HD theory of scalar fields ( see j7], [8], [9], [10], [11] and refer- 
ences therein). It should be mentioned that scalar fields play very important 
role in the theory of inflation of the universe. So, inflation scalar fields can be 
candidates of dark matter of the universe. HD scalar fields are also present 



1 



in SUSY field theories in extra dimensions (see, for example, [12J). Although 
HD field theories lead to ghosts, it was shown that problems with negative 
probabilities and S-matrix unitarity can be solved The purpose of this 
paper is to formulate the forth order equation for scalar fields in the form 
of the first order relativistic wave equation, to obtain the Lagrangian, con- 
served currents, the quantum-mechanical Hamiltonian, and to perform the 
quantization. 

The paper is organized as follows. In Sec. 2, we derive the first order 
relativistic wave equation for scalar fields in the 10-dimensional matrix form, 
the relativistically invariant bilinear form and the Lagrangian formulation. 
In Sec. 3, the canonical energy-momentum tensor and the electromagnetic 
current density are obtained. Dynamical and non-dynamical components of 
the wave function are separated and the quantum-mechanical Hamiltonian 
is found in Sec. 4. In Sec. 5 projection operators extracting solutions of 
field equations for definite energy and different mass states of particles are 
obtained. The canonical quantization of scalar fields with two mass states is 
performed, and the propagators of scalar fields were found in the formalism of 
the first order in Sec. 6. We discuss results obtained in Sec. 7. In Appendix, 
some useful products of matrices are given. The system of units h = c = 1 
is chosen, Greek and Latin letters run 1, 2, 3, 4 and 1, 2, 3, correspondingly. 

2 Scalar Field Equation of the Forth Order 

Consider the forth order field equation describing scalar particles possessing 
two mass states [7], [9]: 



where d 2 = d 2 ., d u = d/dx v = (d/dx m ,d/d(it)). It is obvious that Eq. (1) 
has two solutions corresponding to mass m x and m 2 . 
Let us introduce the 10-dimensional wave function 




(1) 



where <po( x ) — < Pn{ x ) 





(2) 
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mi + ni2 ' 

The function 4>(x) represents the direct sum of two scalars (p(x), (p(x), and 
two four- vectors (f^x), <fn(x). 

Introducing the elements of the entire matrix algebra e A ' B (see, for ex- 
ample [15] ) with matrix elements and products 



(e M > N ) 
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$MA$. 
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where A, B, M, N = 0,0, fi, fi, Eq. (1) can be cast in the form of the first 
order equation 



-0,M 



as 
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where 



a 



m\ + m\ 



rrii + ^2 m\m2 
and there is a summation over all repeated indices. We define 10-dimensional 
matrices as follows: 



r 0,M 



r 0,At 



<T£°" - 



rn 



(6) 



7 10 = e°'° + e°'° + + e^, (7) 

where Jio is a unit 10-dimensional matrix. Then Eq. (5) becomes the rela- 
tivists wave equation of the first order: 



{py.dy. + m) <j>{x) 



0. 



(8) 



So, we reformulated the higher derivative equation for scalar fields (1) in the 
form of the first order Eq. (8). Now, one can apply general methods to 
investigate the first order matrix equation [13] . The spectrum of the particle 
mass of Eq. (8) is given by m/Aj, where Aj are the eigenvalues of the matrix 
p&. It is not difficult to verify that the matrix 



satisfies the matrix equation 

a m? + nin 9 m?m! , . 

^ - / 1 \A + 7 ^^I A = °> 9 

(m 1 + m 2 ) 2 (m 1 +m 2 ) 4 

where 

A = £ 0,0 + £ 0,0 + £ 4,4 + £ 4,4 (1Q) 

is the projection operator extracting four- dimensional subspace, A 2 = A. As 
the eigenvalue of the matrix A is unit, we obtain from Eq. (9) eigenvalues of 
the matrix p 4 as follows: 

A 1 = ±-n* a 2 = ±^^. (id 

mi + m 2 mi + m 2 

So, positive masses of scalar particles described by the first order Eq. (8) are 
ra/Ai = mi, m/A 2 = m 2 . 

The Lorentz group generators in the representation space are 



= e w - + e"' u - (12) 
and obey the commutation relations 

[Jfiv, Jap] — + S^J^a — S u /3J fla — 5^ a J u p. (13) 

The relativistic form-invariance of Eq. (8) follows from the relationship 

[px, J^u] = huPv - huPn- (14) 

One may verify with the help of Eq. (4), (6), (12) that Eq. (14) is valid. 
The Hermitianizing matrix r\ has to obey the relations [33] 

VPm = -ptv + , VP* = ptv + O = !> 2, 3). (15) 

We obtain 



_ £ Q,0 _ £ 0,0 _ a ( m l + m 2) fm,m _ £ 4,4 

m ^ 



(mi + m 2 ) 



£ rn,m _|_ £ m,rn g^'^ 



(16) 



m 

The matrix 77 is the Hermitian matrix, 77 + = rj. Introducing the matrix 
<f)(x) = <f) + (x)r] (<f) + (x) is the Hermitian-conjugate wave function), one obtains 
from Eq. (8) the "conjugate" equation 

0(x) (p^ - m) = 0. (17) 
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The Lagrangian is given by the standard equation 
1 



(x) (p^ + m) (j)(x) - (j){x) [Pfi dp to) <j)(x 



(18) 



so that the relativistically invariant bilinear form is <f>(x)<f>(x) = (x)rj(j)(x) . 
With the aid of Eq. (2), (3), (16), the Lagrangian (18) becomes 



C 



1 



a 
2 



mi + VTL2 I 2TOiTO2 

+ (av) <p 



'f&lp + (ay) <p 



+ TOim 2 V5 V 5 f- 



(19) 



Up to total derivatives, which do not influence on an equation of motion, the 
Lagrangian (19) takes the compact form 



C 



1 



(dj) v ip*) (d^d u tp) 



m 1 m 2 (m 1 + m 2 

+ (ml + to|) (d^if*) (d^ip) + m\m 2 2 Lp*Lp 
One may verify that the Euler-Lagrange equation [T3] 



(20) 



dC 



d, 



dC 



dnd v - 



dC 



dtp* »d(d^*) v d{d»d v v*) 



0. 



(21) 



for the higher derivative Lagrangian (20), leads to the field equation (1). 

3 The Energy- Momentum Tensor and Elec- 
tromagnetic Current 



With the help of the general expression [16 



T, 



fJ,U 



-d u (b{x) + d u 6(x) — t-^= — - - 8 iaJ C, 



(22) 



we obtain from the Lagrangian (18) the canonical energy-momentum tensor 
T r> = 9 ( 9 ^( x )) P^( x ) ~ 7^P( x )p^<P{x)- (23) 



It was taken into account that C = for functions (j){x), <p(x) satisfying Eq. 
(8), (17). Using Eq. (2)-(4), (6), one finds from Eq. (23) the expression as 
follows: 

Tp,v = 1 +m ) {( m i + m s) [<p*d^d u ip - {d^ip*) d v (p] - <f*dfj,d u d 2 <f 

(24) 

- (cPcp j d^cp + d v d\ + (d„<p*) d^d 2 cp\ + c.c, 

where c.c. means the complex conjugate expression. The energy density and 
the momentum density are given by S = T 44 , P m = iT m4 . 
The electric current density is 



- dC dC 

Ux)= '^ (x) j(^))- rnm^ x)] - (25) 

Replacing Eq. (18) into Eq. (25), one obtains the electric current density 

j M (x) = uj)(x)p fl (f)(x), (26) 
so that d^j^x) = 0. With the help of Eq. (2)-(4), (6), we find 

in = r~, A (ml + mf) [{dptp*) <p - <p*d^} 

m-iinoimi + m 2 l v ' 



+ip*d^ip - (fyflV ) if + (d 2 <p*) dtf> - (d^*) d 
For the real scalar fields, ip = cp*, the electric current vanishes 



V 



(27) 



4 Quantum-Mechanical Hamiltonian 

Introducing an interaction of scalar particles with electromagnetic fields by 
the substitution — > = —ieA^, where is the four- vector potential 
of electromagnetic fields (the minimal electromagnetic interaction), Eq. (8) 
may be represented as follows: 



ip 4 d t (f)(x) 



p a D a + m + eA p 4 



cf>(x). (28) 
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Let us introduce two auxiliary functions 

E(x) = A<p{x), Q{x) = U(f)(x), 



(29) 



where the projection operator A is given by Eq. (10), and the projection 
operator IT is 

n = 1 - A = e m ' m + £"*>"*, (30) 
so that E(x) + Q(x) = 4>(x). Acting on Eq. (28) by the operator 

(mi + m 2 ) 2 



2 2 
mfm 2 



"P4 



m{ + m 2 2 - (mi + m 2 ) 2 p\ 



and taking into consideration Eq. (9), we obtain 

id t E(x) = eA E(x) 



(mi + m 2 ) 
H 5 — o — P4 



(31) 



2 2 



m\ + ml- (mi + m 2 ) 2 P4 (p a D a + m) 
Multiplying Eq. (28) by the operator II, one finds 

Ilp a D a E(x) +mtt(x) = 0. 



(32) 



We took into account here that IL94 = 0, Ilp m Il = 0. It follows from Eq. 
(31), (32) that the dynamical components of wave function are given by the 
E(x), and non-dynamical components by the Q(x). The S(x) possesses four 
components, and the auxiliary function Q(x) has six components. To sepa- 
rate the dynamical and non-dynamical components, we express the auxiliary 
function Q(x) from Eq. (32), and replace it into Eq. (31). As a result, one 
obtains the equation in the Hamiltonian form: 



id t Z(x) = HE(x), 



H = eA + 



(mi + m 2 y 
m\m2 



m 1 +m 2 - (mi + m 2 ) p 4 



P4 



x (p a D a + m) 1 Up m D 



m 



(33) 



(34) 



Two components of the function E(x) correspond to the state with the mass 
mi (with positive and negative energy) and the other two - to the state with 
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the mass m 2 . With the help of products of matrices, given in Appendix, the 
Hamiltonian (34) becomes: 



H = eA - ame 4 >° - (m x + m 2 ) e ' 4 + e ' 4 + ml e 4 >° - e 4 >° 



+ (5 4 '° + 5 4 '°)^ 

mi + m 2 



(35) 



We have omitted the linear term in derivatives L = (e A ' m + e 4 ' m ) D m in the 
Hamiltonian, because LE(x) =0. Taking into consideration Eq. (2), (4), (10), (35), 
we can rewrite Eq. (33) in the component form 

id t cp = eA ip - (mi + m 2 ) y? 4 , id t (p = eA ip - (m 1 + m 2 ) (p 4 , 

id t (p 4 = eA tp A -m(p-\ ; D^tf, (36) 

mi + m 2 

id t (pi = eA^ifi — crmip + rrup H D^. 

mi + m 2 

The system of equations (36) may be obtained from Eq. (1), (3), after the 
exclusion of components 

1 Q ~ 1 Q _ 

¥m = ; a m <p, ip m = ■ d m ip. 

m i + m 2 m i + m 2 

According to Eq. (33) only components with time derivatives enter Eq. (36), 
and describe the evolution of fields in time. 



5 Mass Projection Operators 

Consider solutions of Eq. (8) with definite energy and momentum for two 
mass states, r = 1,2, in the form of plane waves: 



ml 



where V is the normalization volume, p 2 = —m 2 . (no summation in index 
r). We imply that four momentum p = (p,ipo) possesses the additional 
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quantum number r = 1, 2 corresponding two masses. The 10- dimensional 
function v T (±p) obeys the equation 

(ip ± m) v T (±p) = 0, (38) 

where p = p^p^. It is convenient to use the normalization conditions 

J ^\x)p^(x)d 3 x = ±1, J $ k \x)p 4 <f>W(x)cPx = 0, (39) 

where <p^\x) = (d)^(x)) rj. Normalization conditions (39) lead to rela- 
tions for functions v T (±p): 

v T (±p) Pll v T (±p) = t^t-, v T (±p)v T (±p) = 1. (40) 

It is not difficult to verify, with the help of Eq. (4), (6) (see Appendix), that 
the minimal equation for the matrix p is 

p - 7 7 \2 P + 7 — 7^2^ = °- ( 41 ) 

(mi + m 2 ) (mi + m 2 ) 

Now we obtain the projection matrices corresponding to definite energy, mo- 
mentum, and quantum number r: 



U T (±p) 



(mi + m 2 ) 4 zp (ip =F m J ^2 



•1 



P + 



(mi + m 2 ) 



(42) 



2mfm2 (m\ — m^mQ, 
The projection operators (42) obey equations as follows: 

(ip ± m) n T (±) = 0, (43) 

n T (± P ) 2 = n r (±p), n r (+p)n T (- P ) = 0, trn T (± P ) = 1. (44) 

Projection matrices (42) can be represented (see jTTj ) as matrix-dyads 

U T (±p)=v T (±p)-v T (±p), (45) 

so that matrix elements of matrix-dyads are 

(v T (±p) ■ v T (±p)) MN = (v T (±p)) M (v T (±p)) N . 

Projection operators (42) extract solutions of Eq. (38) for definite energy 
and different mass states of particles. Eq. (42), (45) allow us to calculate 
matrix elements of different processes of interactions of scalar particles in the 
covariant form. 
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6 Field Quantization 

One can obtain the momenta from Eq. (18): 

dC 



7T [X] 



ir[x) 



x = ^0(^)P4, 

d(d (j){x) 2 



dC 



-Pa<I>{x), 



(46) 



(47) 



where the fields 4>(x), 4>{x) are independent "coordinates". The Hamiltonian 
density is given by the equation 



H = ir(x)d (f)(x) + (d (f)(x)^ n(x) - £ 



= -(t){x)pid Q (t){x) - - (<9 o 0(x)J Pa<I>(x), 

It follows from Eq. (23) that the value Ti = T 44 is the energy density. 
In the quantized theory the field operators can be written as follows 



(48) 



r(z) = E <p^ +) (z) + b r,P<Pr \x) 



(")( 



(49) 



where positive and negative parts of the wave function are given by Eq. (37). 
The creation and annihilation operators of particles a+ , a T)P , and creation 
and annihilation operators of antiparticles b+ p , b TtP obey the commutation 
relations: 

i a T, P , at>, P >] = Srr'Spp', [a T , P , a T ',p'] = [a+ p , a^y] = 0, 

[ & r, P , = 5 TT ,5 PP ,, [6 TiP , ivy] = [b+ p , fc+, p/ ] = 0, (50) 

[a T)P , 6 T / )P /] = [a T)P , &+, p ,] = [a+ p , 6 r /y] = [a+ p , &+, p ,] = 0. 

With the aid of Eq. (48)- (50), and normalization condition (39), we obtain 
the Hamiltonian 



H = Hd 3 x = $> (at,p a r,p + KpKp) ■ 

J T, P 



(51) 
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It is not difficult to find from Eq. (49)- (50) commutation relations as follows: 

[<PtM(x),(J)tn(x')} = [0 rM (x),0 Tjv (x / )] = 0, (52) 

[4>tm(x), 4>tn( x ')} = N rMN (x, x'), (53) 
N rMN (x,x') = N+ MN (x,x') - N~ MN (x,x'), 

Kmn(x^') = E K +) W) M (SV))^ (54) 
One obtains from Eq. (49): 



JV 



2 

N tMn( x , x ') = J2 K(±P)) M {v r {±p)) N exp[±ip(x - x% (55) 

Taking into consideration Eq. (42), (45), and the relation p 2 = —m 2 , we find 
from Eq. (55): 



i(mi + m 2 ) \ m T p (ip =p m) 
2pQVmm 2 m 2 {m\ — m 2 m 2 ) 



p 2 + 



■i 



(mi + m 2 y 



MN 



X exp [±ip{x - X')} = ± 2 2, 4 2 ' ( 56 ) 

mm\m 2 [m% — m\m 2 ) 



x 



4 

77C 



OA 



(m 1 +m 2 ) 2 \)mn p 2p V 

With the help of the singular functions p2] 

A +(» = E 77-77 exp(zpx), A_(x) = E 77-77 exp(-ipx), 
p ^PoV p ^PoV 

A (x)=i{A + {x)-A„{x)), 
we obtain from Eq. (54), (56) 

(mi + m 2 ) 4 m 2 T (ppd^) (p M <9 M - m) 



N tMN (x,x') = -i 



mm 2 m 2 (m^ — m 2 m 2 ) 



(57) 
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X 



\ Ao(x-x'). 

) MN 



(mi + m 2 ) 2 

For the points x and x', which are separated by the space-like interval (x — 
x') > 0, the commutator [4)m(x), 4> N (x')} equals zero due to the properties 
of the function A (a;) [16J. If t = t', [0 M (x, 0), ^(x', 0)] = N rMN (x - x', 0), 
and the function N rMN (x — x', 0) can be obtained from Eq. (57) with the 
help of equations 

d 2 n A (x)\ t=0 = 0, ^A (a;)| t= o = 0, d A (x)\ t=0 = 5(x), (58) 

where n = 1,2,3,.... One may verify, using Eq. (41), the validity of the 
equation 

(^ (1 + m)iV T ± (x, a :')=0. (59) 

The propagator of scalar fields (the vacuum expectation of chronological 
pairing of operators) can be defined in our formalism as 



(T(j) TM {x)(f) TN (y))o = K MN {x - y) 

= 0(x - y ) N+ MN (x -y)+9(y - x ) N~ MN (x - y), 
where 9{x) is the well known theta-function. We obtain from Eq. (56): 

(mi + m 2 ) 4 m 2 . (p^<9 M ) (p^ - m) 



(60) 



(T(f) TM (x)(j) TN (y)}o 



(61) 



x 



ml 



MN 



A c (x-y), 



(mi + m 2 ) 2 
and the function A c (x — y) is given by 

A c (x - y) = 9 (x - y ) A + (x -y) + 9(y Q - x ) A_(x - y). 



(62) 



Propagators (61) are finite only for mi 7^ m 2 . It is seen from Eq. (61) that 
propagators have different signs for r = 1 and r = 2. This means that one 
of states of the scalar field is the ghost [9]. 
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7 Conclusion 



We have formulated the scalar field equation with higher derivatives in the 
form of the 10-component first order relativistic wave equation. This equa- 
tion describes scalar particles possessing two mass states. It should be noted 
that the second order equation for scalar fields with one mass state is formu- 
lated in the 5-component matrix form [18], [19]. The relativistically invari- 
ant bilinear form, and the Lagrangian were obtained, and this allowed us to 
find the canonical energy-momentum tensor and density of the electromag- 
netic current by the standard procedure. We found the quantum-mechanical 
Hamiltonian by the separation of dynamical and non-dynamical components 
of the wave function. The wave function entering the Hamiltonian equation 
possesses four components to describe two mass states of scalar fields with 
positive and negative energies. The first order relativistic wave equation as 
well as the Hamiltonian equation are convenient for different applications. 
The density matrix (matrix-dyad) found can be used for calculations of elec- 
tromagnetic possesses. The first order formalism allowed us to quantize HD 
scalar fields in a simple manner. The HD scalar field theory considered can 
by applied for a model of inflation of the universe, where one of states of 
the scalar field is identified with Quintessence and another, the ghost, with 
Phantom [U] . 
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With the help of Eq. (4), we find useful products of matrices entering the 
Hamiltonian (34): 



m 



p\p m 



III 



mi + VTL2 



^^A,m i _4,m _4,m 

ae ' + e — e 



,mi + 77T.2 

p 4 p m Ilp n = 5, 



a z - 1 ) e 4 ' m + ae 4 ' m - ae 4 ' m - e 4 ' m 



m 



p\p m Ilp n = 5 n 



mi + 777,2 



mi + m 2 

3 



£ 4,0 _ £ 4,0 _ a£ 4,0 



£ 4,0 +(7£ 4 J 0_ (T£ 4 ) + ^_ (T 2N £ 4 1 

p 4 n = o. 



(63) 

(64) 

(65) 

(66) 
(67) 
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Using the definition p = p^p^, we obtain: 



p 2 = ^ ( ae°>° + - £ °> Cl 
mi + m 2 

m 



p^le^-e^-ae^ 



(68) 



A o mp 

V = ; 

mi + m2 



mi + m 2 



{mi + ^ 2 ) 2 

, mp 2 cr o m 2 p 2 

p = p 

mi + m 2 {nil + m 2 ) 2 



P 2 



(69) 



(70) 



f = K + m2) ^V - mV 2 p, (71) 
(mi + m 2 ) (mi + m 2 ) 
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